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Abstract 

Lagrangians linear in velocities were analyzed using the fractional 
calculus and the Euler-Lagrange equations were derived. Two examples 
were investigated in details, the explicit solutions of Euler-Lagrange equa- 
tions were obtained and the recovery of the classical results was discussed. 

PACS: ll.10.Ef. Lagrangian and Hamiltonian approach 

1 Introduction 

The mathematical idea of fractional derivatives, which goes back to the seven- 
teenth century, has represented the subject of interest for various mathematicians 
[1, 2, 3]. The fractional calculus, which means the calculus of derivatives and inte- 
grals of any arbitrary real or complex order is gaining a considerable importance in 
various branches of science, engineering and finance [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. 

As it was observed during the past three decades, the fractional differential 
calculus describes more accurately the physical systems [4, 5, 6, 12]. Many applica- 
tions of fractional calculus amount to replacing the time derivative in an evolution 
equation with a derivative of fractional order. This is not merely a phenomeno- 
logical procedure providing an additional fit parameter. One of the problems 
encountered in this field is what kind of fractional derivatives will replace the 
integer derivative for a given problem [1, 2, 3, 13]. Depending on the specified 
physical situation different authors have applied different derivatives [6, 9]. 

Nonconservative Lagrangian and Hamiltonian mechanics were investigated by 
Riewe within fractional calculus [14, 15]. Besides, Lagrangian and Hamiltonian 
fractional sequential mechanics, the models with symmetric fractional derivative 
were studied in [16, 17] and the properties of fractional differential forms were 
introduced [18, 19]. 

Recently, an extension of the simplest fractional variational problem and the 
fractional variational problem of Lagrange was obtained [20]. A natural and 
interesting generalization of Agrawal's approach [20] is to apply the fractional 
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calculus to the constrained systems [21, 22]. In the present work, we apply the 
concept of fractional calculus to the Lagrangian with linear velocities. 

The plan of this paper is as follows: 
In Sec. 2 the Riemann-Liouville (RL) fractional derivatives were briefly presented. 
In Sec. 3 the fractional Lagrangians with linear velocities were constructed and 
their corresponded Euler-Lagrange were analyzed. Sec. 4 was devoted to our 
conclusions. 



2 Riemann - Liouville fractional derivatives 

In this section the definitions of the right and left RL derivatives as well as their 
basic properties are briefly presented. 

RL fractional derivatives [3, 20] are defined as follows 



W"-f(^y(-s) /< T - t)"'°' 1 f(T)dT, (2) 

where the order a fulfills n — 1 < a < n and T represents the Euler's gamma 
function. The first derivative is called the RL left fractional derivative and in (2) 
the expression of the RL right fractional is presented. If a is any integer, the 
relations to the usual derivative are obtained as follows 

Under the assumptions that f(t) is continuous and p > q > , the most 
general property of RL fractional derivatives can be written as 

a n p t ( a npf(t)) = a nr q f(t). (4) 

For p > and t > a we obtain 

JD?( a D t -7(*)) (5) 

which means that the RL fractional differentiation operator is a left inverse to 
the RL fractional integration operator of the same order. The relation (5) is 
called the fundamental property of the RL fractional derivative. In addition, the 
fractional derivative of a constant is not zero and the RL fractional derivative of 
the power function {t — a) v is given by 
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where v > — 1. The normal derivatives and a D\ commute only if f^\a) = 
0, j — 0, 1, . . . , n — 1 is fulfilled and two RL fractional derivative operators Df 
and D ( 5 commute only if 



7(*) 







(7) 



and 



t=a 







(8) 



The above properties of RL fractional derivatives lead us to the conclusion that 
there are many substantial differences from the usual derivatives and therefore the 
solutions of the differential fractional equations contain more information than 
the classical ones. 



3 Fractional Euler- Lagrange equations for 
Lagrangians with linear velocities 

Let Jlq 1 , ■ • ■ , q n \ be a functional of the form 

b 

j L(t,q\---,q\ a D t V, • • • , a D«g", t D? q\ • • • , t Dfc») dt (9) 

a 

defined on the set of functions q l (t), i = l,---,n which have continuous left 
RL fractional derivative of order a and right RL fractional derivative of order 
j3 in [a, b] and satisfy the boundary conditions q l (a) = q l a and q l {b) = q\. A 
necessary condition for Jlq 1 , ■••,<?"] to admit an extremum for given functions 
q l (t), i — 1, • ■ • ,n is that q l (t) satisfy Euler-Lagrange equations [20] 

In the following we consider the Lagrangian with linear velocities 

L = a s (<t)q i -V(j), (11) 

where a,j (q l ) and V {q l ) are functions of their arguments. 

The fist step is to construct the corresponding fractional generalization of 
the Lagrangian given by (11). The fractional Lagrangian is not unique, in other 
words there are several possibilities to replace the time derivative with fractional 
ones. The requirement is to obtain the same Lagrangian expression if the order 
a is 1. Having in mind the above considerations, for < a < 1, we propose two 
fractional Lagrangians. The first one is as follows 
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L> = a 3 (q*) a VW-V(q*) . (12) 
From (10) and (12), the corresponding Euler-Lagrange equations emerge as 

^JW+.D^-^O. (13) 

The second Lagrangian is given by 

U = -a, (g*) t T>«q> - V (j) . (14) 
Using (10) and (14) the corresponding Euler-Lagrange equations become 

^ W - + »D- aiM + ^ = 0. (15, 
3.1 Examples 

A. To illustrate our approach, let us consider the following Lagrangian 

L = q l q 2 - q 2 q l - (q 1 - q 2 )q 3 (16) 

which is a gauge invariant [23]. In this case we proposed the corresponding 
fractional Lagrangian to be as 

U = („D t y ) q 2 - ( a B?q 2 ) q l - (q l - q 2 )q 3 . (17) 
Using (13), the Euler-Lagrange equations corresponding to (17) become 

q^q 2 , - a V?q 2 -q 3 + t V«q 2 = 0, J^g 1 + q 3 - t D^ = 0. (18) 
The solution of (18) is given as follows 

q 1 = q 2 , (19) 

q 3 = (- a r>? + t r>«)q 1 . (20) 

From (19) and (20) we conclude that the classical solution is obtained if a — > 1. 
B. Let us consider the second Lagrangian given by 

L = q\ 2 + qY-V(q 2 ,q 3 ,q*), (21) 

where V(q 2 ,q 3 ,q 4 ) = ^[(g 4 ) 2 ~~ 2g 2 g 3 ]. We observe that (21) is a second class 
constrained system in Dirac's classification [21]. 
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We propose the fractional generalization of (21) to be as follows 

U = -[( t D?gV + GDgtfV + V(l 2 ' <?i 9 4 )] • (22) 
From (15) and (22) the Euler-Lagrange equations are given by 

aVIq 2 = , (23) 

i D^ 1 + g 3 = 0, (24) 

a Btq A + q 2 = 0, (25) 

t D^ 3 - g 4 = . (26) 
From (23), we conclude that the solution for q 2 (t) has the form 

q 2 (t) = C 1 (t-a) a - 1 . (27) 

From (25) and (27) an equation for q 4 (t) is obtain as follows 

aD^^-C^-o) - 1 . (28) 

The solution of (28) has the form 



q \t) = C 2 {t - a)"- 1 -Jjl-J(-T + tr-^r - ar-'dr . (29) 



t 



Using (26) and (29), the solution of q 3 (t) becomes 



q 3 (t) = Csib-tr-' + ^-Ji-t + rr-'i-a + rr-'dr 



Ci 

2 



t 

b T 



r(a) 



y (-* + r)"- 1 J(t- 77) Q_1 (r7 - a^drjdr . (30) 



Finally, the equation (24) together with (30) give the solution for g 1 (t) as 
follows 



c, " 



«'(*) = ftfi-ir 1 -^ /(-( + T)«- 1 (6-T)<>- 1 dT 

1 (a) J 



/(-( + t)"- 1 I (-t + ^""'(-o + nr^drjdr 

t T 
b b V 

J (a- t)"- 1 J {a- n)"- 1 j(r- v )°-\r, - a^dTA/d^l) 

\ / t cr a 



T( a f 
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Here C±, C 2 , C 3 and C 4 are constants. If a— » 1, a — ► 0, 6 — ► 1, then the standard 
solutions 

q\t) = ^-^f+C'j+C'^ q \t) = C' A , g 3 (t) = f t 2 -C 3 i+c4 g 4 (i) = -C' 4 t+C' 3 

(32) 

are recovered if we redefine the constants from (27), (29), (30) and (31) as 

C 1 — C4 — C3 ^ + "3"' C*2 = — 2^ + C*2 + C*3) C3 — C2, C A = C\. 

4 Conclusion 

The Lagrangians with linear velocities were investigated using left and right RL 
fractional derivatives. The corresponding fractional Lagrangians were proposed 
and the fractional Euler-Lagrange equations were obtained. Although the frac- 
tional Lagrangians contain only left or right derivatives, both derivatives are in- 
volved in Euler- Lagrange equations and both played an important role in finding 
the solutions. The exact solutions of the Euler-Lagrange equations were obtained 
for two examples corresponding to first and second-class constrained systems. 
A "gauge invariance" was reported for the first example. The solutions of the 
investigated examples depend on the limits a and b and the limiting procedure 
recovered the standard results. 

5 Acknowledgments 

This work is partially supported by the Scientific and Technical Research Council 
of Turkey. One of the authors (D. B.) would like to thank R. Hilfer, O. Agrawal 
and N. Engheta for providing him important references and M. Naber and M. 
Henneaux for interesting discussions. 

References 

[1] K. B. OLDHAM K. B. and SPANIER J., The Fractional Calculus, (Academic 
Press, New York) 1974. 

[2] MILLER K.S.and ROSS B., An introduction to the Fractional Integrals and 
Derivatives-Theory and Applications, (Gordon and Breach, Longhorne, PA) 
1993. 

[3] PODLUBNY L, Fractional Differential Equations, (Academic Press) 1999. 

[4] MAINARDI F., Fractional calculus:Some basic problems in continuum and 
statistical mechanics, in A.Carpinteri, F. Mainardi (Eds.), Fractals and Frac- 
tional Calculus in Continuum Mechanics, ( Springer- Verlag, New York) 1997. 



6 



[5] MAINARDI F. and PAGNINI G. , Appl. Math. Comp., vol. 141, no. 1 

(2003) 51. 

[6] HILFER R. (Ed.) , Applications of Fractional Calculus in Physics, (Word 
Scientific Publishing Company, Singapore, New Jersey, London and Hong 
Kong) 2000. 

[7] RIDA S. Z. , AMA El-Sayed, Apl. Math. Com. 147 (2004) 29. 

[8] ROSU H.C., MADUENO L. and SOCORRO J., J. Phys. A 36 (2003) 1087. 

[9] AGRAWAL P. O. , Nonlinear Dynamics 29 (2002) 145. 

[10] ENGHETA N., Micr.Opt.Tech.Lett. vol. 22, no.4 (1999) 237. 

[11] ANH V. V. and NGUYEN C. N., Finance Stoch. 5, no.l (2001) 83. 

[12] BLUTZER R. L. and TORVIK P. J., J. Rheology 30 (1996) 133. 

[13] HILFER R., Chem. Phys. 284, no.(l-2) (2002) 399. 

[14] RIEWE F., Phys. Rev. E53 (1996) 1890. 

[15] RIEWE F., Phys. Rev. E55 (1997) 3582. 

[16] KLIMEK M., Czech. Journ. Phys., 52, no.ll (2002) 1247. 

[17] KLIMEK M., Czech. Journ. Phys., 51, no.12 (2001) 1348. 

[18] SHEPHERD C. and NABER M., J. Math. PhysA2, no. 5 (2001) 2203. 

[19] SHEPHERD C. and NABER M., math-ph/0301016. 

[20] AGRAWAL O. P., J. Math. Anal. Appl. 272 (2002) 368, 

AGRAWAL O. P., Trans. ASME, J. Math. Anal. Appl, 68, no. 2 (2001) 
339. 

[21] HANSON A., REGGE T. and TEITELBOIM C, Constrained Hamiltonian 
systems (Academia Nationale dei Lincei, Rome) 1976; 

SUNDERMEYER K., Constrained Dynamics, Lecture Notes in 
Physics, vol. 169 ( Springer- Verlag, New- York) 1982); 

HENNEAUX M. and TEITELBOIM C, Quantization of Gauge Sys- 
tems, (Princeton University Press) 1992. 

[22] BALEANU D. and GULER, Y., Nuovo Cimento B, 115, no.3, (2000) 319. 

[23] SUGANO R. and KIMURA T., Phys. Rev. D41 (1990) 1247. 



7 



